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Abstract 

The Invaded Cluster (IC) dynamics introduced by Machta et al. [Phys. 
Rev. Lett. 75 2792 (1995)] is extended to the fully frustrated Ising model 
on a square lattice. The properties of the dynamics which exhibits numerical 
evidence of self-organized criticality are studied. The fluctuations in the IC 
dynamics are shown to be intrinsic of the algorithm and the fluctuation- 
dissipation theorem is no more valid. The relaxation time is found very short 
and does not present critical size dependence. 
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I. INTRODUCTION 



Recently a new type of cluster Monte Carlo (MC) dynamics, the Invaded Cluster (IC) 
dynamics, based on invaded percolation has been introduced by Machta et al. for the ferro- 
magnetic Ising model. |l| The IC dynamics is based on the Kastelein-Fortuin and Coniglio- 
Klein (KF-CK) cluster formulation of the Ising model [Q] and has been shown to be even 
more efficient than the Swendsen-Wang (SW) dynamics in equilibrating the system at 
the critical temperature. The IC dynamics has the advantage that the value of the critical 
temperature does not need to be known a priori. In fact the dynamics itself drives the 
system to the critical region as in self-organized critical (SOC) systems. 

The aim of this paper is to extend the IC dynamics to frustrated systems where KF-CK 
clusters percolate at a temperature Tp, higher than the critical temperature T^. In particular 
we will consider the fully frustrated (FF) Ising model on a square lattice where it has been 
shown numerically that ksTp/J ^ 1.69 (A;^ is the Boltzmann constant and J is the 
strength of the interaction) and = 0. We will use two definitions of clusters. The first 
definition is based on the KF-CK clusters. In this case the IC dynamics leads to a SOC 
percolating state at temperature Tp corresponding to the percolation of KF-CK clusters. 
The second definition is more general ^ and reduces to the cluster of Kandel, Ben-Av and 
Domany (KBD) |^ in the FF Ising model. In this case the IC dynamics leads to a SOC 
state at the thermodynamical critical temperature Tc = 0. 

In sec. II we review the definition and the results of IC dynamics on ferromagnetic Ising 
model. Then we extend the IC dynamics to the FF Ising model using the KF-CK cluster in 
sec. Ill and the KBD clusters in sec. IV where we study also the equilibrium relaxation of 
the proposed dynamics. 
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II. IC DYNAMICS FOR FERROMAGNETIC ISING MODEL 



The rules which define the IC dynamics for the ferromagnetic Ising model are very 
simple. Let us start from a given spin configuration. As first step all the pairs of nearest 
neighbour (nn) spins are ordered randomly. Then, following the random order, a bond is 
activated between the nn spin pairs only if the two spins satisfy the interaction. The set 
of activated bonds partitions the lattice into clusters which are classified every time a new 
bond is activated. When one of the clusters spans the system the procedure is stopped 
and the spins belonging to each cluster are reversed all together with probability 1/2. The 
spin configuration obtained can be used as starting point for the next application of the 
dynamical rule. Successively, the algorithm has been generalized to different stopping rules 
|0J^ and to different model (the Potts model and the Widom-Rowlinson fluid 0). 

The rational which is behind this rule is the well known mapping between the ferro- 
magnetic Ising model and the correlated bond percolation. p| In this framework bonds are 
introduced in the ferromagnetic Ising model between parallel nn pairs of spins with prob- 
ability p = 1 — exp(— 2/9J) where P = l/^ksT). The clusters (KF-CK clusters), defined 
as maximal set of connected bonds, represent sets of correlated spins and percolate exactly 
at the critical temperature, Tp = T^. The well known SW cluster dynamics uses these 
clusters to sample very efficiently the phase space at any temperature. In this dynamics at 
each MC step the KF-CK clusters are constructed and the spins belonging to each cluster 
are reversed all together with probability 1/2. This produces a new spin configuration on 
which a new cluster configuration can be build up. The IC dynamics is very similar to 
SW dynamics and differs only on the way clusters are constructed. Since the dynamic rule 
introduced by Machta et al. builds up clusters which, by definition, percolate through the 
system, it is expected that the average properties of the clusters are the same as the KF-CK 
clusters at the critical temperature. In fact, in Ref. [0] it has been shown that the ratio of 
activated bond to satisfied interactions, in the large L limit {L is the linear lattice size), is 
very close to the critical probability Pc = 1 ~ exp(— 2J/(A;bTc)) with which KF-CK clusters 
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are constructed and it has a vanishing standard deviation. Furthermore, the estimated mean 
energy also tends in the large L limit to a value very close to the critical equilibrium value 
E{Tc) with the finite size behaviour expected in the ferromagnetic Ising model at the critical 
point. On the other hand the energy fiuctuation C = {E"^) — {E)"^ is not related to the spe- 
cific heath. In fact, it has been found that C diverges linearly with L and not logarithmically 
as the specific heat does. The latter result points to the fact that the IC dynamics does 
not sample the canonical ensemble in finite volume and, therefore, energy fiuctuations and 
specific heat are not any longer related by the fiuctuation-dissipation theorem. Since the 
fiuctuation C^^^/L^ for L i— > oo, it is generally assumed that in the thermodynamic 
limit the IC dynamics is equivalent to the canonical ensemble even if a rigorous proof is still 
lacking. 

III. IC DYNAMICS FOR FF ISING MODEL WITH KF-CK CLUSTERS 

We now extend the IC dynamics to frustrated systems where the KF-CK clusters per- 
colate at a temperature Tp which is higher than the critical one T^. In particular we will 
consider the FF Ising model on square lattice where Tp ~ 1.69 ^ and Tc = 0. The FF 
Ising model is defined by the Hamiltonian 

H = -JJ2{e,,S,S,-l) (1) 

(ij) 

where Si takes the values ±1 and assumes the value —1 on even columns and is +1 
otherwise. 

In the simplest extension of the IC dynamics we introduce bonds at random between 
spins satisfying the interaction (i.e. eijSiSj = 1) until a spanning cluster is found. As in the 
ferromagnetic case, at this point the spins belonging to each cluster are reversed all together 
with probability 1/2. The procedure is iterated until equilibrium is reached. We have done 
simulations performing measurements over 5 ■ 10^ MC sweeps after discarding the first 10^ 
for equilibration for systems with sizes L ranging from 16 to 300, and over 3.75 ■ 10^ MC 
sweeps after discarding the first 750 for L = 350. 
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In Fig.l.a we show as functions of the system size L the ratio of activated bonds A*";, to 
satisfied interaction A^^ = — E /2 (where E is the energy), (/) = Nh/Ng, and in Fig.l.b 
the density of energy e = E/L"^. Assuming the size dependences (/) = (/)oo ~ A/L and 
e = e^ + A/ log(L) we find (/)oo = 0.698 ± 0.008 and eoo = 1-24 ± 0.01, which are very close 
to the values Pc = 1 - e^'^/'^^ 0.694 and e(Tp) = 1.234 at the temperature Tp ^ 1.69. The 
errors on (/) and e decrease with increasing L, denoting that the (/) and the e distributions 
become sharp when L i-^ oo . 

The results obtained show that the dynamics has not driven the system into the critical 
thermodynamical state at T = 0, but to the percolation critical state at T = Tp. In 
order to study the convergence as function of L we have extracted for the IC dynamics at 
each L an effective temperature Tjc{L) (see Tab.l) by using (/) = 1 — exp(— 2/T7C'), then 
we have compared the IC energy density to the analytical energy density e{Tjc) at Tic 
(Fig.l.b). This analysis clearly show that the L i-^ oo limit is reached in the IC dynamics 
very slowly. In particular the energy is systematically larger than that obtained by SW 
dynamics. Furthermore, since the dynamics build up, by definition, percolating clusters the 
mean cluster size S/ of IC clusters diverges in the following way = AL>/'^f where 

Ip/^p = 1.78 ±0.11 is a good approximation of the exponent "^p/vp = 1.792 of Random Bond 



Percolation |jTT| and A = 0.022 ± 0.013 (Fig. 2). This result is in excellent agreement with 
the behaviour found for KF-CK clusters for both the critical exponent and the prefactor 
A. ^ This critical behaviour drives also energy and magnetization fluctuations in a critical 
regime, i.e., (E^) - (E)'^ ~ L^-^ and (M^) - (M)^ ~ L^-^ (Fig.3). This result is in strong 
contrast with the behaviour of specific heat and magnetic susceptibility at T = Tp which 
behave as L^. In fact Tp is not a critical thermodynamic temperature and both specific heat 
and magnetic susceptibility are finite at T = Tp. This again stresses the fact that the IC 
dynamics does not sample the canonical ensemble in finite volume and shows that the spin 
configurations visited by the IC dynamics with a finite probability correspond to a much 
larger range of energies and magnetization than in any ordinary dynamics. 

It is interesting to note that the straightforward application of the IC dynamics to more 



complex systems as the Ising spin glass (SG) model could allow to sample very efficiently the 
equilibrium spin configurations of the system at T = Tp where the KF-CK clusters percolate 



(in 2d SG Tp ^ 1.8 §, in 3d SG Tp ~ 3.95 ||T2[) with a dynamics which exhibits SOC. 
However since Tp is usually much larger than the spin glass critical temperature Tsg (in 2d 
TsG = 0, in 3d Tsg — l-H 0), any other dynamics does not suffer for slowing down. 



IV. IC DYNAMICS FOR FF ISING MODEL WITH KBD CLUSTERS 

We ask now how to build up an IC dynamics able to drive a system to the frustrated 
thermodynamical critical temperature. From the analysis and the argument given for the 
ferromagnetic case we understand that the cluster definition needs to be modified in such 
a way that i) the frustrated system can still be mapped on the corresponding correlated 
percolation model and ii) the clusters percolate at the thermodynamical critical temperature. 
The first conditions is always satisfied by the KF-CK clusters |T5[ but not the second one. 



A general procedure to construct such clusters in a systematic way was suggested in Ref . . 
In particular for the FF model this procedure leads to the cluster algorithm proposed 
by Kandel, Ben-Av and Domany (KBD). |[16| To define the clusters in KBD dynamics 
one partitions the square lattice in a checkerboard way and chooses randomly one of the 
two patterns (Fig. 4a). For each plaquette if three of four spin pairs are satisfied (a single 
plaquette in the FF model can have either one or three satisfied spin pairs) one activates 
bonds between the two spin pairs satisfying the interaction and facing each other (Fig. 4b) 



with a probability p = 1 — exp(— 4/T). Numerically it was showed [0 that the clusters 
obtained with such procedure percolates at the thermodynamical critical temperature T^ = 
0, with critical exponents z/p = 1 and 7p = 2 to be compared with the thermodynamical 
critical exponents u = 1 and 7 = 3/2. 

With this idea in mind we propose the following invaded cluster dynamics. From the 
checkerboard partition we order randomly all the square plaquettes belonging to the chosen 
pattern. Then plaquettes are tested in this order to see how many spin pairs are satisfied and, 
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in plaquettes with three spin pairs satisfying the interaction, we activate bonds between the 
two spin pairs satisfying the interaction and facing each other. Every time a pair of bonds 
is activated the cluster structure change and the occurring of a spanning cluster is checked. 
As in the previous cases when the first cluster percolates the cluster evolution is stopped 
and a new spin configuration is obtained by reversing the spin belonging to each cluster all 
together with probability 1/2. The dynamics proposed is related to the algorithm introduced 
by KBD exactly as the IC dynamics by Machta et al. is related to the SW algorithm. 

We have tested the proposed algorithm performing measurements over 10^ MC sweeps on 
square lattices of sizes L ranging from 16 to 200 and over 5-10'^ MC sweeps on L = 250-i-400, 
after discarding the first 10^ for equilibration. 

In Fig. 5. a we show the results of our simulations which give a situation very similar to that 
obtained in the ferromagnetic case. In fact, the ratio of plaquettes with activated bonds Ni, to 
plaquettes with three satisfied interactions Ns = (3L^— i?)/4, (/) = Nb/Ng, reaches the value 
Pc = 1 for large system size L. At the same time the effective temperature Tjc (see Tab. 2), 
obtained by (/) = l—exp{—4:/Tjc), converges rapidly to the limit value = 0. In Fig. 5. a we 
shows alsopp = 1— exp(— 4/Tp(L)) where Tp{L) is the estimated |]T^ percolation temperature 
of KBD clusters in a square FF system of size L. The mean energy per spin (Fig.S.b) e tends, 
within the numerical precision obtained, towards the thermodynamic value at the critical 
temperature e(T = 0) = 1 with the expected size dependence e — e^c = AL~^ + BL^^ . The 
fit of the data gives eoo = 1.016 ±0.002. The errors on (/) and e go to zero with increasing L. 
The estimated mean cluster size exponents '^pjvp = 1.2 ± 0.2 (Fig. 6) does not coincide with 
the expected values 7p/z/p = 2. [|1^ We explain this result with the very slow convergence as 
function of L of the percolation quantities obtained with the KBD-clusters (see for example 



Fig. 8 of Ref. [0). Thus we expect to recover the right behavior only for sizes fairly large. 
We have studied the energy and magnetization fluctuations also in this case (Fig. 7) and 
found that {E^) - {Ef ~ L^-^ and (M^) - {Mf ~ L^-^. These exponents, as in the previous 
cases do not agree with those expected for specific heat and magnetic susceptibility. As in 
the ferromagnetic case we obtain that energy and magnetization fluctuations are larger in 



the IC dynamics than in the canonical ensemble. 

We have also studied the equilibrium relaxation of the magnetization of the proposed IC 
dynamics 

^ {M{trM{t'+tr)-{M{t')r ... 

{M{t'Y) - {M{t')Y 

where the time t is measured in MC steps. As shown in Fig.8 (j){t) vanishes in few MC steps. 
The integrated autocorrelation time r defined as 

1 tm 

r-^ + Y^m (3) 

is reported in Tab. 3 for different values of L. The dependence on 1^, is extremely week since 
T reaches a plateau very quickly. The value obtained r ~ 1.6 is lower than that obtained in 
the KBD dynamics tkbd ^ 2.4 and show a weak tendency to decrease with increasing L. 



V. CONCLUSIONS 

We have shown how the IC dynamics introduced by Machta et al. for the ferromagnetic 
Ising model can be extended to the FF Ising model on a square lattice. The straightforward 
extension with KF-CK clusters shows that the IC dynamics leads to a Self Organized Critical 
(SOC) percolating state at the percolation temperature Tp and can be used to produce 
equilibrium spin configurations at temperature different from the thermodynamical critical 
temperature, actually at the percolation temperature Tp. The dynamics is characterized 
by intrinsic diverging fluctuations and the fluctuation-dissipation theorem is no more valid. 
The extension with KBD clusters, whose percolation point coincide in the large L limit 
with the critical point of the FF system, has properties very similar to those obtained in 
the ferromagnetic model: it drives the system to the critical region without a previous 
knowledge of the critical temperature and gives a reasonably good estimation of the average 
energy at the critical point. The estimated integrated autocorrelation time is smaller than 
that obtained in the KBD dynamics. We have also stressed that the extension has been 
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possible since there exists a percolation model into which the FF square Ising model can be 
exactly mapped. The extension to other frustrated systems, such as spin glass, in principle 
can be done using the systematic procedure suggested in Ref. 0. The computation have 
been done on DECstation 3000/500 with Alpha processor. 
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TABLES 







TABLE I. 


KF-CK 


clusters: 


numerical (^sduiat 


cs of Tj( 


(L).m 






L 


16 


50 


64 


80 


100 


150 


200 


250 


300 


350 


Tic{L) 


1.98 


1.77 


1.75 


1.73 


1.72 


1.703 


1.695 


1.689 


1.686 


1.682 


error 


0.04 


0.02 


0.01 


0.01 


0.01 


0.007 


0.007 


0.005 


0.004 


0.004 


TABLE IL KBD clusters: numerical estimates of Tic{L).[lQ\ 


L 


16 


50 


64 


80 


100 


150 


200 


250 


300 


400 


Tic{L) 


1.05 


0.76 


0.71 


0.67 


0.64 


0.58 


0.54 


0.51 


0.49 


0.46 


error 


0.15 


0.17 


0.12 


0.12 


0.14 


0.16 


0.19 


0.19 


0.25 


0.46 



TABLE in. Integrated autocorrelation time r for IC dynamics with KBD clusters and tkbd 
for KBD dynamics 



L 


100 


200 


300 


400 


T 


1.74 ± 0.02 


1.73 ± 0.02 


1.67 ±0.02 


1.60 ± 0.02 


TKBD 


2.800 ± 0.008 


2.48 ±0.02 


2.56 ±0.01 


2.37 ± 0.03 
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FIGURES 

FIG. 1. KF-CK clusters: ratio of activated bonds to satisfied interaction (/) (a) and energy 

density e (b) versus system size L for a square FF system; squares are the IC dynamics output and 
the superimpose fines are the plot of the fit functions (/) = {f)oo — A/ L (with (/)oo = 0.698±0.008 
and A = 1.2±0.2) and e = eoo+^/log(L) (with = 1.24±0.01 and A = 0.62±0.07) for L > 100; 
circles are the energy density values of a standard MC dynamics at the estimated temperature Tjc 
in Tab.l; dashed line is the analytical energy density for the model at Tjc; dotted line is the 
asymptotic value at T = 1.69 ~ T^; arrows show the asymptotic values estimated by the fits; 
where not shown the errors are included in the symbols; all quantities are dimcnsionless.[10] 

FIG. 2. KF-CK clusters: mean cluster size per spin S/L'^ for IC dynamics for a square FF 
system with L = 16-=- 350. We show the fit over data points for L > 150 with the best fit parameters 
given in the text {A = 0.022ib0.013 and jp/vp = 1.78±0.11); the errors are included in the symbols; 
all quantities are dimensionless.[10] 

FIG. 3. KF-CK clusters: density of energy fluctuation (a) and of magnetization fluctuation 
(b) for a square FF system; squares are the fluctuations for IC dynamics; circles are the fluctuation 
values of a standard MC dynamics at the estimated temperature Tjc in Tab.l; dotted line is the 
analytical values for the model at Tjc] the estimated behaviour are {E'^) — {E)'^ /L? ~ L^''^ and 
(M^) — (M)^/L^ ~ L^-^ (the error on the exponents is on the last given digit); where not shown 
the errors are included in the symbols; all quantities are dimensionless.[10] 

FIG. 4. Square FF lattice: the solid (dashed) lines are ferromagnetic (antiferromagnetic) inter- 
actions, (a) Checkerboard partition of a square lattice: the plaquettes with (without) a dot give a 
pattern; (b) an example of plaquette with three satisfled interactions: in this case we activate the 
two vertical bonds. 
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FIG. 5. KBD clusters: ratio of plaquettes with activated bonds to plaquettes with 3 satis- 
fied interactions (/) (a) and energy density e (b) versus system size L for a square FF system; 
squares are the IC dynamics output and the superimpose hues are the plot of the fit functions 
(/) = (/)^ - A/L - B/L'^ (with (/)oo = 1.000 ± 0.003, A = 0.1 ± 1.2 and S = 9 ± 100) and 
e = eoo-A/L-B/L^ (with eoo = 1.016±0.002, A = 4.3±0.5 and B = 1.0±0.7) for L > 100; stars 
are the vahics oi Pp = 1 — exp(— 4/Tp(L)) with Tp{L) percolation temperature of KBD clusters in a 
square FF system of size L (from Ref.[17]); circles are the energy density values of KBD dynamics 
at the estimated temperature Tjc in Tab. 2 (the errors are asymmetric because they are derived 
from the indetermination on Tjc)', dashed line is the analytical energy density for the model at 
Tjc; dotted line is the asymptotic value at Tc = 0; arrow shows the asymptotic value estimated by 
the fit; where not shown the errors are included in the symbols; all quantities are dimensionless.[10] 

FIG. 6. KBD clusters: mean cluster size per spin S/L'^ for IC dynamics for a square FF 
system with L = 16 400. We show the linear fit ln(S'/L^) = A + 'Jp/i^p In(-L) over data points for 
L > 250 with A = —0.8 ± 1.5 and jp/vp = 1.2 it 0.2; the errors are included in the symbols; all 
quantities are dimensionless.[10] 

FIG. 7. KBD clusters: the density of energy fluctuation (a) and of magnetization fluctu- 
ation (b) for a square FF system; squares are the fluctuations for IC dynamics; circles are the 
fluctuation values of a standard MC dynamics at the estimated temperature Tjc in Tab. 2 (the 
error are asymmetric because are derived from the indetermination on Tjc); dashed line is the 
analytical values for the model at Tjc', the estimated behaviour are {E'^) — {E)"^ / ~ L*^'^ and 
(M^) — (M)^/L^ ^ L}-^ (the error on the exponents is on the last given digit); where not shown 
the errors are included in the symbols; all quantities are dimensionless.[10] 

FIG. 8. Magnetization correlation functions in a square FF system for the KBD dynamics (for 

L = 100 at T ~ Tp{L) ^ 0.45 and for L = 400 at T = Tp ~ 0.342 [17]) and for IC dynamics with 
KBD clusters (for L = 100 and L = 400). Time t is measured in MC steps. 
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